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We present results for diflerent observables in weak decays of pseudoscalar and vector mesons 
with a heavy c or fe quark. The calculations are done in a nonrelativistic constituent quark model 
improved at some instances by heavy quark effective theory constraints. We determine pseudoscalar 
and vector meson decay constants that within a few per cent satisfy fvMv/ fpMp = f, a result 
expected in heavy quark symmetry when the heavy quark masses tend to infinity. We also analyze 
the semileptonic B ^ D and B — » D* decays for which we evaluate the different form factors. 
Here we impose heavy quark effective theory constraints among form factors that are not satis- 
fied by a direct quark model calculation. The value of the form factors at zero recoil allows us 
to determine, by comparison with experimental data, the value of the \Vcb\ Cabbibo-Kobayashi- 
Maskawa matrix element. From the B —> D semileptonic decay we get \Vcb\ = 0.040 ± 0.006 in 
perfect agreement with our previous determination based on the study of the semileptonic A;, — > Ac 
decay and also in excellent agreement with a recent experimental determination by the DELPHI 
Collaboration. We further make use of the partial conservation of axial current hypothesis to de- 
termine the strong coupling constants gs* BniO) — 60.5 ±1.1 and go* Dir(O) — 22.1 ± 0.4. The ratio 

= ( gB*Bn{0) .fs'VMD )/( go'DniO) fo'VM^ ) = 1.105 ± 0.005 agrees with the heavy quark 
symmetry prediction of 1. 

PACS numbers: 12.15.Hh,12.39.Hg,12.39. Jh,13.20.Fc,13.20.He 

I. INTRODUCTION 

In systems with a heavy quark with mass much larger than the QCD scale (Aqcd) a new symmetry, known as heavy 
quark symmetry (HQS) [ESS Hi arises. In that limit the dynamics of the light quark degrees of freedom becomes 
independent of the heavy quark flavor and spin. This is similar to what happens in atomic physics where the electron 
properties are approximately independent of the spin and mass of the nucleus (for a fixed nuclear charge). HQS can 
be cast into the language of an effective theory (HQET)0 that allows a systematic, order by order, evaluation of 
corrections to the infinity mass limit in inverse powers of the heavy quark masses. HQS and HQET have proved very 
useful tools to understand bottom and charm physics and they have been extensively used to describe the dynamics 
of systems containing a heavy c or 6 quark For instance, all lattice QCD simulations rely on HQS to describe 

bottom systems 0. 

In a recent publication we have studied the A° h^l^Di and S+Z^i// reactions in a nonrelativistic 

quark model. The detailed analysis of the different form factors showed how a direct nonrelativistic calculation does 
not meet HQET constraints and we had to improve our model imposing HQET relations among form factors. Our 
calculation allowed for a determination of the \Vcb\ Cabibbo-Kobayashi-Maskawa (CKM) matrix element given by 
\Vcb\ = 0.040 ± 0.005lg'QQ2 in good agreement with a recent determination by the DELPHI Collaboration \Vcb\ = 
0.0414 ± 0.0012 ± 0.0021 ± 0.0018 ^10] . What we intend to do here is a study of different weak observables of 
pseudoscalar and vector mesons with a heavy c or 6 quark in a non relativistic quark model improved at some points 
with HQET constraints. Weak observables are of great interest as they help to probe the quark structure of hadrons 
and provide information to measure the CKM matrix elements. 

In the case of mesons with a heavy quark HQS leads to many model independent predictions. For instance in the 
HQS limit the masses of the lowest lying (s-wave) pseudoscalar and vector mesons with a heavy quark are degenerate. 
Nonrelativistic quark models satisfy this constraint: the reduced mass of the system is just the mass of the light quark 
and the spin-spin terms, that distinguishes vector from pseudoscalar, are zero if the mass of the heavy quark goes 
to infinity. HQS also predicts that the masses and leptonic decay constants of pseudoscalars and vector mesons are 
related via fp Mp — fv My, relation that is also satisfied in the quark model in the HQS limit. If one looks now at 
the form factors for the semileptonic B ^ D and B D* decays HQS predicts relations among different form factors 
that are also met by the quark model in the HQS limit. The question is to what extent the deviations from the HQS 
limit evaluated in the quark model agree with the constraints deduced from HQET. In addition we will make use of 
these HQET constraints to improve the quark model results and thus come up with reliable predictions. 

The paper is organized as follows: in section^ we introduce the meson wave functions and interquark potentials 
we shall use in this work. In section HTll we analyze the leptonic decays of pseudoscalar and vector B and D mesons, 
determining the different decay constants. In section llVl we study the form factors for the semileptonic B ^ Dlv and 
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B D* Iv decays. In section^we evaluate the strong coupling constants qe'E-k and go'D-K- Finally in section IVTI 
we end with the conclusions. The paper also contains three appendices where we collect the expressions for the matrix 
elements that are needed for the evaluation of different observables. 

Apart from lattice QCD and QCD sum rules (QCDSR) calculations with which we shall compare our results and that 
will be quoted in the following, the different observables analyzed in this work have been studied in the quark model 
starting with the pioneering work of Ref. within a non relativistic version, to continue with different versions of 
the relativistic quark model applied to the determination of decay constants [l^ IT^ IT^ ITp. Il^ IitI ll^ [l9l form 
factors and differential decay widths 19, 21, 22, 23,^4,2^,261,23, Isgur-Wise functionsIlllillillalM^ ,33i ,3| 
or strong couphng constants [s^, IIB iS IS ^ • 



II. WAVE FUNCTION AND INTERQUARK INTERACTIONS 



For a meson M we use the following expression for the wave function 



0.1,0.2 



q, ai pi = 



q, a.2 P2 



V2 



-P+p 



(1) 



where P stands for the meson three momentum and A represents the spin projection in the meson center of mass, ai 
and a2 represent the quantum numbers of spin (s), flavor (f) and color (c) 



a = (s,/,c) 



(2) 



of the quark and the antiquark, while Ef-^, pi and Ef^, p2 are their respective energies and three-momenta, m/ is the 
mass of the quark or antiquark with flavor /. The factor (—1)^"*^ is included in order that the antiquark spin states 
have the correct relative phase^ . The normalization of the quark and antiquark states is 

{a' p'\ap) ^S^,, „ i2nf 2E Sip ' - p) 



(3) 



Furthermore, (f>aTa'2 ^P) momentum space wave function for the relative motion of the quark-antiquark system. 

Its normalization is given by 



and, thus, the normalization of our meson states is 



NR 



M, X' P'\M,XP 



NR 



6x'.x{2tt)H{P'-p) 



(4) 



(5) 



For the particular case of ground state pseudoscalar (P) and vector (V) mesons we can assume the orbital angular 
momentum to be zero and then we will have 



1 



}.(P) 



V3 

^'5,,,,, H) ^^^f.^{\p\) YooQ) (1/2, 1/2,0; Si, S2,0) 
(-1) "^ooQ) (1/2, 1/2,1; Si, S2, A) 



(6) 



^ The list of references is by no means exhaustive. 

^ Note that under charge conjugation (C) quark and antiquark creation operators are related viaC Ca{p) = (— dL{p). This means 
that the antiquark states with the correct spin relative phase are not d]y.{p) |0) = | <?, 01 p) but are instead given by (— d]^(p) |0) = 
(-1)5- |g, cp). 
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where (ji, J2, Js; "t-i, "22, ma) is a Clebsch-Gordan coefficient, Yqq = l/y/in is tlie / = m = splierical harmonic, and 
^^f^^f^ilpl) the Fourier transform of the radiai coordinate space wave fimction. The phases are introduced for later 
convenience. 

To evaluate the coordinate space wave function we shall use several interquark potentials, one suggested by Bhaduri 
and collaborators [s^ (BHAD), and four suggested by Silvestre-Brac and Semay |40|,i4i| (AL1,AL2,AP1,AP2). The 
general structure of those potentials in the quark-antiquark sector is 

J r niinij rr^ SmiiTij V ^ tt^Xq j 

with a the spin Pauli matrices, rrii the constituent quark masses and 

XQ{mi,mj)=A\ ■ — (8) 

\ mj + nij J 

The potentials considered differ in the form factors used for the hyperfine terms, the power of the confining term [p — 1, 
as suggested by lattice QCD calculations ^2|> or p = 2/3 which gives the correct asymptotic Regge trajectories for 
mesons '43^ ) , or the use of a form factor in the one gluon exchange Coulomb potential. All free parameters in the 
potentials have been adjusted to reproduce the light (tt, p, K , K*, etc.) and heavy-light {D, D*, B, B* , etc.) meson 
spectra. They also lead to precise predictions for the charmed and bottom baryon (Ac^f,, Sc.b, 5]*^, Sc,b, j^, S* j^, 
flc^b and ri* ^) masses liol, and for the semileptonic A° k^l'vi and E'^l~Pi decays. 

We will use the above mentioned interquark interactions to evaluate the different observables. This will provide us 
with a spread of results that we will consider, and quote, as a theoretical error to the averaged value that will quote 
as our central result. 



III. LEPTONIC DECAY OF PSEUDOSCALAR AND VECTOR B AND D MESONS 

In this section we will consider the purely leptonic decay of pseudo scalars {B, D) and vector {B* , D*) mesons. 
The charged weak current operator for a specific pair of quark flavors /i and /2 reads 



(Ci, Si), (C2, S2) 



*ai(0)7p(l- 75)^02(0) 



(9) 



with ^'q.j a quark field of a definite spin, flavor and color. The hadronic matrix elements involved in the processes 
can be parametrized in terms of a unique pseudoscalar fp or vector fv decay constant as 



P,P) 



J 



fih 

A 1^1 



(0) 



p,p 



V,XP) = {0\ (0) V,\P)= e^-^^ ( P )Mvh 



'V ^ 



where the meson states are normalized such that 



M, \' P' \ M,\Pj = Sy^x (27r)3 2E S{P' - p) 



(10) 



(11) 



In the first of Eas. (|10|) is the four-momentum of the meson, while in the second My and eji'^'' ( P ) are the mass and 
the polarization vector of the vector meson. In both cases /i and /2 are the flavors of the quark and the antiquark 
that make up the meson. 

Concerns about the experimental determination of the pseudoscalars decay contants have been raised in Ref. [45j . 
There the effect of radiative decays was analyzed concluding that for B mesons the decay constant determination 
could be greatly affected by radiative corrections. In the vector sector, and as rightly pointed out in Ref. 0|, the 
vector decay constants are not relevant from a phenomenological point of view since B* and D* will decay through 
the electromagnetic and/or strong interaction. They are nevertheless interesting as a mean to test HQS relations. 

For mesons at rest we will obtain 



fv 



Mv 



(0| J, 



P,0 

y, 00 



(12) 
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with Mp the mass of the pseudoscalar meson. In our model, and due to the different normalization of our meson 
states, we shall evaluate the decay constants as 



fp = ~i 



2 



2 



P,0 



NR. 



y. 00 ) 

/ NR 



(13) 



The corresponding matrix elements are given in appendix IXI 

The results that we obtain for the different decay constants appear in Tables and ^] Starting with /^i and /^^ 
our results are larger that the ones obtained in the lattice by the UKQCD Collaboration [H^l or the ones evaluated 
using QCD spectral sum rules (QSSR) j^. Not only the independent values are larger but also the ratio fn Hp i s 
larger in our case. On the other hand our results are in better agreement with other lattice determinations [Sfl l5(i| . 
They also compare very well with the experimental measurements of /d and jn, in Refs. lisL 11^ Isollsil l53ll b eing 
our I Jd ratio in very good agreement with the value obtained using recent CLEO Collaboration data "Bill. As 
for /b and fs^ , we find a very good agreement in the case of fs^ between our results and the ones obtained in the 
lattice or with the use of QSSR. For fs our result is smaller and then also our ratio fBs/fs is larger. 





fo [MeV] 


fo, [MeV] 


foJfD 


This work 


243li7 


34111 


i-4ilS:^^ 


Experimental data 








CLEO 
ALEPH \m\ 
OPAL [50] 
BEATRICE [51j 
E653 
BES [53] 


202 ± 41 ± 17 [42 280 ± 19 ± 28 ± 34 [4^ 

— 285 ± 19 ± 40 

— 286 ± 44 ± 41 

— 323 ± 44 ± 12 ± 34 

— 194 ± 35 ± 20 ± 14 
3711111 ±25 — 


— 
— 


Lattice data 








UKQCD [54] 
Fermilab Lattice [5^ (Preliminary) 
M. Wingate et al. [56] 


2251^3 ±21 


229(3)1?^ 

263I9 ± 24 
290 ± 20 ± 29 ± 29 ± 6 


1-11(1)11 


QCD Spectral Sum Rules 








S. Narison [51] 


203 ± 23 


235 ± 24 


1.15 ±0.04 




/fl [MeV] 


/ss [MeV] 


fsJfB 


This work 




239I7 


i.54l°:«^ 


Lattice data 








UKQCD [54] 
M. Wingate et al. [56] 
Lattice world averages 


200 ± 30 [^ 


220(6)lig 
260 ± 7 ± 26 ± 8 ± 5 
230 ±30 [^ 


1.13(1)11 
1.16 ±0.04 [^ 


QCD Spectral Sum Rules 








S. Narison [^ 


207 ± 21 


240 ± 24 


1.16 ±0.04 



TABLE I: Pseudoscalar fp decay constants for B and D mesons 
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This work 
UKQCD [54] 


9 1+0-9 

8.6(3)1^ 


fi 5+°-^ 

0.0_o.4 

8.3(2)15 


1 .-,+0.06 
J^-*^-0.05 

1.04(1)1^ 




/s* [MeV] 


/sj [MeV] 




This work 
UKQCD [54] 


35.6111 
28(l)li 


23.0+1° 
25(ljl^ 


1.10(2)1^ 



TABLE II: fv = Mv / fv for B* and D* mesons 



For the vector meson decay constants wc obtain the values 

fa- = 223li^MeV" fo* = 326+^^ 



Ib' = 151 



h23 - 
-19^ 

+ 15 
-13 



Jb, = 236lliMey 



(14) 



which are very much the same as the values obtained for the decay constants of their pseudoscalar counterparts. This 
almost equality of pseudoscalar and vector decay constants is expected in HQS in the limit where the heavy quark 
masses go to infinity where one would have Q 



fvMv = fpMf 



My = Mp 



(15) 



Our decay constants satisfy the above relation within 2%. On the other hand UKQCD lattice data show deviations 
as large as 20% for D mesons [53| . 

In order to compare the values of the vector decay constants with lattice data from Ref. [E^l we give in Table lll| 
the quantity fv — My/ fv- We find good agreement for //j. and /s; but not so much for the other two. Also our 

ratios fo'lfoi and /s* //s* are larger than the ones favored by lattice calculations. 
On the other hand the ratio 



1.138: 



0.008 



is in very good agreement with the expectation in Ref. [G^I where they would get 1.05 ~ 1.20 for that ratio. 



(16) 



IV. SEMILEPTONIC DECAY OF B INTO Dlv AND D*lv 

In this case the strong matrix elements are parametrized as 

D; P' I J^''(0) \B,P) ( D, P' I J{,\{Q) I B, P 



^/MbMd 



D*, XP' I J^''(0) I B, P 



iv + v')^ h+iw) + {v-v')^ h^{w) 



^ fiiyaf: 



,p (e(^)''(P')) v^v'^hv{w) 
+ 1 (ef^) {P'))* -v {v^ hA, (w) + v'^ hA, (w)) 



(17) 



(18) 



where v = P/Mb and v' ~ P' /Mjj^ are the four velocities of the initial B and final D, D* mesons, w = v - v' ^ and 
e^ivap is the fully antisymmetric tensor with £0123 = +1. 



^ ui is related to the four momentum transferred square via = A4g + — 2 w Mg Mjj^ £>• 
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w 


P+ 


P- 


Pv 


Pa, 


Pa^ 


Pa, 


1.0 


2.6 


-5.4 


11.9 


-1.5 


-11.0 


2.2 


1.1 


-0.3 


-5.4 


8.9 


-3.8 


-10.3 


-0.2 


i.Z 


—3.1 


—5.3 


6.1 


—5.9 


—9.8 


— 2.5 


1.3 


-5.6 


-5.3 


3.5 


-7.9 


-9.3 


-4.6 


1.4 


-8.0 


-5.2 


1.1 


-9.7 


-8.8 


-6.6 


1.5 


-10.2 


-5.2 


-1.1 


-11.5 


-8.4 


-8.5 


1.59 


-12.1 


-5.1 











TABLE III: QCD corrections in % as evaluated in Ref. ^ 



In the limit of infinite heavy quark masses mc^rnt — > oo HQS reduces the six form factors to a unique universal 
function ^{w) known as the Isgur-Wise function U 

h+{w) = hA,{w) = hAsiw) = hv{w) = £,{w) (19) 
h.{w) ^ hA^iw) =0 (20) 

Vector current conservation in the equal mass case implies the normalization 

= 1 (21) 
Away from the heavy quark limit those relations are modified by QCD corrections so that one has 

h, (w) = {a, + P, (w) + 7, (w) + 0{l/ml,)) ^w) (22) 
The aj are constants fixed by the behavior of the form factor in the heavy quark limit 

a+ = aAi = 01A3 — cty — 1 

Q_ ^aA2^0 (23) 

The different Pj account for perturbative radiative corrections |6l| while the 7^ are non perturbative in nature and 
are proportional to the inverse of the heavy quark masses j63 |. At zero recoil (w = 1) Luke's theorem [6^ imposes 
the restriction. 



7+(l)=7Ai(l) = (24) 

so that power corrections to /i+(l) and /i^j(l) are of order 0{l/ml ^). In Tables HTn and HVI we collect the values for 
the different f3j and 7^ in the full interval of w values allowed in the two decays. These two tables have been taken 
from Ref. .62] . 



w 


7+ 


7- 


7v 


lAi 


7A2 


7A3 


1.0 


0.0 


-4.1 


19.1 


0.0 


-23.1 


-4.1 


1.1 


2.7 


-4.1 


20.7 


2.9 


-21.4 


-0.7 


1.2 


6.2 


-4.1 


23.1 


6.5 


-19.8 


3.4 


1.3 


10.5 


-4.2 


26.3 


10.7 


-18.3 


8.0 


1.4 


15.3 


-4.4 


30.0 


15.4 


-17.0 


13.0 


1.5 


20.6 


-4.5 


34.3 


20.5 


-15.8 


18.5 


1.59 


25.7 


-4.7 











TABLE IV: Power corrections 'jj{w) in % as evaluated in Ref. [6^ 
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A. B —> Dlu decay 



Let us start with the B ^ Dlv case. In the center of mass of the B meson and taking P' — —q — —\q \ k in the z 
direction we will have for the form factors ft.+ (w) and h^{w)* 

h+iw) = \ fy"(|<f I) + (EniWl) - Md)) 

y/2MB2MD \ \q\ J 

h^{w) = \ (v'm + (EniWl) + Mn)) (25) 

V2Mb2Md V \<l\ ) 

where EB{\q\) = ^JlvFj^T^ and V^{\q\) (/i = 0, 3) is given by 

V^{\q\) = ( -\q\ k I (J^'nO) I B, O) (26) 
In our model V^^(|9 |) is evaluated as 

V^^m = ^2MB2EDi\q\) ^^^A-kli|(J^')^(0)| i?, 0^ 



NR 



(27) 



which expression is given in appendix FBI 

In the case of equal masses mf, = uic vector current conservation demands that 



/i+(l) = 1 ; h^{w)=0 (28) 

In this limit we find that ^+(1) = 1 so that our value for /i+(l) complies with vector current conservation. On the 
other hand h-(w) ^ violating vector current conservation. 

In the upper left panel of Fig. ^ we show the values of hj^iw) and h-{w) for the B ^ D transition as obtained from 
Eqs. H25I I27|l with the use of the ALl interquark potential. The values for h-{w) are not reliable. Actual calculation 
shows that they are of the same size as the deviations from zero that one computes in the equal mass case. To improve 
on this what we shall do instead is to use the form factor h+{w) and Eq. 122|) to extract ^{w) (we shall call it ^+{1^1) ) 
and from there we can re-evaluate h-{w) with the use of Eq. (|22l) . The results appear in the upper right panel of 
Fig. n] where we also show the lattice results for ^ (w) obtained by the UKQCD Collaboration in Ref . [6J| . We find 
good agreement with lattice data. Finally in the lower panel of Fig. ^ we show the different £,+ {w) obtained with the 
use of the different interquark potentials. As we see from the figure all ^+(w) are very much the same in the whole 
interval for w. 

The slope at the origin of our Isgur-Wise function is given by 



^2 _ 1 d^+iw) 



0.35 ± 0.02 (29) 



£,+ iw) dw 

small compared to the lattice value of = O.Slt.W extracted from a best fit to data. 

1. Differential decay width 

Neglecting lepton masses the differential decay width for the process B ^ DID is given by 

dT Gl 



dw 487r3 



\Vcb\^ Ml (i«2 - lf/^{MB + MdY Fl{w) (30) 



where Gp = 1.16637(1) x 10 ^ GeV ^ |73 is the Fermi decay constant, Vcb is the CKM matrix element for the b ^ c 
weak transition, and Fd{w) is given by 



Foiw) = 



1 - r 

h+M - ——h_{w) 
1 + r 



(31) 



In this case to is related to |(j | via |<? | = Mjj \/ w"^ — 1 
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FIG. 1: Upper left panel: h+{w) and h-{w) for the B ^ D transition as obtained from Eqs. 12511271 using the ALl interquark 
potential. Upper right panel: h+ {w) as before, ^+ (w) obtained from the values of h+ (w) using Eq. 1221 , h- [w) obtained from 
using Eq. 1221 . We also show the UKQCD lattice data by Bowler et al. |64|. Lower panel: Different £,+{w) obtained 
with the above procedure for the different interquark interactions. Lattice data is also shown. 



with r Mo/Mb- 

In Fig. 121 we show our calculation for Fd{w) \Vcb\ obtained with the ALl interquark potential and using three 
different values of \Vcb\ corresponding to the central and extreme values of the range for \Vcb\ favored by the Particle 
Data Group (PDG), \Vcb\ = 0.039 0.044 jT^I- We also show the experimental data for the decays B~ ^ D^lv and 
5° £)+ I V obtained by the CLEO Collaboration [g^, a fit to CLEO data using the form factors of Boyd et al. [6^ . 
and the experimental data for the decay D'^lv obtained by the BELLE Collaboration [g^- Our results are 

larger than experimental data for w > 1.2. Our total integrated width will thus be larger than the experimental one 
for any reasonable value of \Vcb\- From our data we extract the slope at w = 1 given by 



1 dFoiw) 



F]j{w) dw 



= 0.38 ±0.02 



(32) 



which is small compared to the values extracted from the experimental data: pjj = 0.76 ± 0.16 ± 0.08 |6y| and 
pjj = 0.69 ± 0.14 IHg obtained from a linear fit to the data, or pjj = 1.30 ± 0.27 ± 0.14 ^ and pjj = 1.16 ± 0.25 |H| 
obtained using the form factors of Boyd et al. Thus, only our results close to w = 1 seem to be reliable. We 

can use our prediction for -Fb(l) to extract the value of \Vcb\ from the experimental determination of the quantity 
I K:fc|i^D(l). Different values of that quantity appear in Table IVl 

Our result for Fd(1) is given by (we do not show the theoretical error which is or the order of lO^**) 

Foil) = 1.04 (33) 

which is in good agreement with other calculations Fd{1) = 0.98 ± 0.07 69], Fd{1) = 1.04 or 
Fnil) = 1.069 ±0.008 ± 0.002 ± 0.025 [z^- From our value for Fd(1) and the experimental values for |Vcb|fb(l) 
we can obtain \Vcb\ in the range 



0.040 ±0.006 



(34) 
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FIG. 2: FD{r,w) \Vcb\ obtained with the ALl interquark potential. Solid line: our results using |Vc6| = 0.0415. Dashed 
line: our results using \Vcb\ ~ 0.044. Dotted line: our results using \Vcb\ = 0.039. Circles: experimental data by the CLEO 
Collaboration ^G^. Dashed-dotted line: fit to CLEO data using the form factors of Boyd et al. 67']. Diamonds: experimental 
data by the BELLE Collaboration Q. 





\V^b\FD{l) 


CLEO Collaboration [66] 


0.0416 


± 0.0047 ± 


0.0037 


BELLE Collaboration [68J 


0.0411 


± 0.0044 ± 


0.0052 



TABLE V: \Vcb\ Fn{l) values obtained by different experiments. 



This result agrees with our recent determination based on the analysis of the Af, — > Ac Z vi reaction from where we got 
IKhl = 0.040 ±0.005 0. 



B. B ^ D* I u decay 



Working again in the center of mass of the B meson and taking P' = —q = — \q\k in the z direction we will have 
for the form factors hv{w), hAi{w), hA2{w) and hA^iw) the expressions 



hv{w) 
hAi{w) 
hA2 {w) 



Mb 



V2 



w + l ^MbMd' 

( Ai%{\q\) En.iW\)4:liW\) 
Mb [ \q\ |gP 



V2Md' 



^/MbMd' 



with ED*{\q\) = v/m|7+1?F, and where v!^*li\q\) and A^*>^{\q\) are given by 

v!::lm) = (D-\X -\q\k\{J^''),iO)\B,6 



(35) 



(36) 
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FIG. 3: Left panel: hviw), hAi{w), hA2{w) and hA^iw) form factors obtained using Eq. I35II . Right panel: Ri{w) and R2{w) 
ratios. In both panels the ALl interquark potential has been used. 



In our model Vl*l{\q\) and A 



(*) ( 



\q\) are evaluated as 



NR 



NR 



D*,\ -\q\k\{.r/)^{())\B, o) 

I NR 



(37) 



with expressions given in appendix ICl 

In the left panel of Fig. |21we show our results for the /iai (w), ft-Aj (w^) and Ha-Xw) form factors, obtained 

with the ALl interquark potential and the use of Eq. In the right panel of the same figure we show the ratios 



hv{w) 

hA^jw) +rhA2iw) 
hAi{w) 



(38) 



where now r = Md* /Mb- These ratios are expected to vary very weakly with w. We find indeed that this is so in 
our case being our values of Ri{w) and R2{w) within 4% of unity. In Table IVll we give now our results for i?i(l) 
and i?2(l) and compare them to different experimental^ and theoretical determinations. We find discrepancies of the 
order of 15 - 33 % for i?i(l) and 13 - 46 % for i?2(l). 

One can understand these discrepancies by evaluating the different ^{w) functions obtained from the form factors 
with the use of Eq. (|22|l and the /3 and 7 coefficients of Neubert given in Tables Iml and im The results appear in the 
upper left panel of Fig.^] One can infer from the figure that our results for hA2 (w) are not reliable. Also we somehow 
miss the correct normalization for hy{l). On the other hand the values of ^Ai{w) and ^Aaiw) are equal within 4% 
and in reasonable agreement with lattice data from Ref. . 

To improve the nonrelatistic quark model prediction, and similarly to what we did in subsection II V Al we will take 
^Aiiw) as our model determination of the Isgur-Wise function ^{w) and we will reevaluate the form factors with the 
use of Eq. 12211 . What we obtain is now depicted in the upper right panel of Fig. 0| In the lower panel we give the 
different ^^j^ (w) obtained with the different interquark potentials. They do not show any significant difference. 

The slope of the ^Ai {w) function at the origin is given by 



= 0.55 ±0.02 



(39) 



^ The experimental results by the CLEO and BABAR Collaborations have been obtained with the assumption that Ri(w) and R.2{w) 
are constants. 
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i?2(l) 


This work 


1.01 ± 0.02 


1.04 ± 0.01 


CLEO [71] 


1.18 ±0.30 ±0.12 


0.71 ±0.22 ±0.07 


Di\rii\ix (^r^reiiminaryj \iZ\ 


l.oZo ± U.UOO ± U.UZO ± U.U^O 


A OOn _l_ A n/1 /I -LA AO A _L A Al Q 

U.yzU ± U.U44 dr U.U^U ± U.Uio 


Caprini et al. 


1.27 


0.80 


Grinstein et al. 


1.25 


0.81 


Close et al. [281 


1.15 


0.91 



TABLE VI: Ri{l) and i?2(l) 



to be compared to the lattice result — O.OSlgg [6J|. In this case we are within lattice errors, but one can not be 
very conclusive due to the large value of the latter in this case. 




0,5 



-- 4v(") 












.-. ^(») 




■ Lattice ^ (m) 





1,2 



CD 



1,4 





FIG. 4: Upper left panel: different £,{w) functions obtained from the hviw), hAiiw), hA2{'w) and liA^iw) form factors using 
Eq. 12211 and the ALl interquark potential. Lattice data by K. C. Bowler et al. from Ref. 16411 are also shown. Upper right 
panel: form factors obtained from 5ai(u)) with the use of Eq. 12211 . Lower panel: Different S,Ai(w) obtained with the different 
interquark potentials. 



Finally in Fig. Owe give the ratio S^^{w)/S^Ai {w) evaluated with the ALl interquark potential. We see the differences 
between the two Isgur-Wise functions are at the level of 3-7%. 
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1,4 



1,2 



0,8 



0,6 



_ ^^((0)/^^ (CO) 



ALl 



1,2 



1,4 



CO 



FIG. 5: Ratio of our two Isgur-Wise functions calculated with the ALl interquark potential. 



1. Differential decay width 
Neglecting lepton masses the differential decay width for the process B — > D* I v is given by [tE 



dw 



487r3 

where F^,. (w) is defined as 



\Ychf (Mb - MD' f Ml. .Jiw-' - 1) [w + 1)' 



Aw 1 — 2wr + \ 



1 



{l-rf 



Fl.{w) 



\ 1-1- 4w l — 2wr+r^ 

The Hj{'w) are helicity form factors given in terms of the and R2{w) ratios as 

w — I 



(40) 



(41) 



Hoiw) = 1 



1 -r 



[l-i?2H] 



Vl — 2wr + 
1 - r 



, w — 1 „ , , 



(42) 



Similarly to Fig. |21 in Fig. Elwe show our results for the quantity F]j.{w) \Vcb\ evaluated with the ALl interquark 
potential and using the values of \Vcb\ corresponding to the central and extreme values of the range for \Vcb\ favored 
by the PDG. We also show the experimental data by the CLEO Collaboration [tJ] for the B~ — > D*°lD reaction 
(squares), and for the — > D*^lv reaction (circles) together with a best fit, and the experimental data by the 
BELLE Collaboration 76] for the i?" D*^ I v reaction (diamonds). We find good agreement with CLEO data for 
small w values. Disagreement starts already at around w = 1.1 where our resuls start to go above the experimental 
data. BELLE data are systematically below our results. 

Also our slope at the origin 



dFu- {w) 



Fd' (w) dw 



= 0.31 ±0.02 



(43) 



is smaller than the value obtained by the BELLE Collaboration pj^, — 0.81 ± 0.12 [t^ using a linear fit to their data. 
All this means that our total width would be larger than the experimental one for any reasonable value of Vcb- On 
the other hand experimentalists are able to extract the value of \Vcb\ F]j<-{1). Different experimental results for that 
quantity appear in Table IVTll 
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> 

* 



0,05 



0,04 



0,03 



0,02 



0,01 



ALl 




o 



This work using |V^^I=0.044 
This work using IV^^I=0.0415 
This work using IV^|^I=0.039 
Exp. D** 1 V (CLEO Coll.) 
Exp. D*" 1 V (CLEO Coll.) 
Fit to CLEO data 
Exp. D*" 1 V (BELLE Coll.) 



1,2 



CO 



1,4 



FIG. 6: FD*{r,w) \Vcb\ obtained with the ALl interquark potential. Solid line: our results using jVcb! = 0.0415. Dashed line: 
our results using IVcbl^ 0.044. Dotted line: our results using \Vcb\ ~ 0.039. Circles and squares: experimental data by the 
CLEO Collaboration [71|. Dashed-dotted line: fit to CLEO Collaboration data. Diamonds: experimental data by the BELLE 
Collaboration iTal. 





\V,t\FD*{l) 


(CLEO Coll.) 


[71] 


0.0431 


± 0.0013 ± 


0.0018 


(DELPHI CoU.) 


fio] 


0.0392 


± 0.0018 ± 


0.0023 


(BELLE CoU.) 


[76] 


0.0354 


± 0.0019 ± 


0.0018 


(BABAR Coll) 


[77] 


0.0355 


± 0.0003 ± 


0.0016 



TABLE VII: \Vcb \ -Fd*(1) values obtained by different experiments. 



Our result for Fo.(l) is given by 

F£,.(l) = /lAi(l) 0.983 ±0.001 (44) 

Comparison with the experimental data for \Vcb \ -Fi5*(l) allows us to extract values for \Vcb\ in the range 

\Vcb\ = 0.0333 - 0.0461 (45) 

One can not be more conclusive due to the dispersion in the experimental data for |Vch| i^D'(l)- From DELPHI data 
alone we would obtain \Vcb\ = 0.040 ± 0.003 in perfect agreement with our determination using the B ^ D reaction 
data. We should also say that our value for Fd.(1) is larger than the lattice determination Fd.(1) = O.OlOjlQ ogg by 
S. Hashimoto et al. normally used by experimentalists to extract their \Vi,b\ values. A new unquenched lattice 
determination of this quantity by the Fermilab Lattice Collaboration is in progress [s^. 



V. STRONG COUPLING CONSTANTS gn-H^ 

In this section we will evaluate the strong coupling constants gH'H-n where H stands for a B or D meson. To this 
end we shall make use of the partial conservation of the axial current hypothesis (PCAC) which relates the divergence 
of the axial current to the pion field as 

d''Jil{x)=Uml<^,-{x) (46) 

where /^r = 130.7 ± 0.1 ± 0.36 MeV [t^ is the pion decay constant, rriTr = 139.57MeV [t^ is the pion mass, and 
$7r- {x) is the charged pion field that destroys a tt~ and creates a 7r+. Using the LSZ reduction formula one can relate 
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the matrix element of the divergence of the axial current to the pion emission amplitude as 

{H, P'|g^'4;(0)|i/*,AP) = -z/.3^^(,\)^jP',P) (47) 

q 

where q = P — P' and A'^I j^^{P' , P) is the pion emission amplitude for the process H* Hn given by^ 

a'^^IhAp'^p) ^-m^HAi') (9^4'^^)) (48) 

The matrix element on the left hand side of Eq. (|^7|i has a pion pole contribution that can be easily evaluated to be 

2 

{H, P'\q^Jil{0)\H*,XP),,on-poie = ~iU J_ , A^h'^,h^{P',P) (49) 

q 

so that we can extract a non-pole contribution 

{H, P'\q^Jil{0)\H\\P)„,„_^,i, = A^^IhAP'^P) 

= -ifn9H'HAq') (q'^s^^Hp)) (50) 

which is the one we shall evaluate within the quark model. For the matrix element on the left hand side of Eq. 150() 
we can use a parametrization similar to the one used in Ea. (|18|l 

{H, P'|g^4^(0)|F*,AP)„„„_p„,, = q^|-»e|,^)(P) {w + 1) hA,H 

+i (e^^) {P ) • «') hA, (w) + hA, {w)) I ^MhMh' (51) 

with the result that 



gH-HTriq ) = ^|(^ + 1)/Mi(w) + W yj^^A2iw) - (w) 



+ 



^ hA, {w) - hA, {w)^ I VMhMh^ (52) 



M 



The evaluation of the form factors is done is a similar way as the one described in subsection IIVBI The results that 
we get for q^ = are 

5D*D.(0) = 22.1 ±0.4 , 5B*B.(0) -60.5 ±1.1 (53) 

to be compared to the experimental determination go* Dtt{^%) — 17.9 ± 0.3 ± 1.9 by the CLEO Collaboration [soj . 
the lattice results go'D-n {ml) = 18.8±2.3l2:S 113 and gB'Biri^) = 47± 5 ± 8 [s^, or a recent determination using 
QCDSR for which gD-DAml) = 14.0 ± 1.5 and gB-B^{'^) = 42.5 ± 2.6 "SS^. Older QCDSR results give smaller 
values for both coupling constants. For instance, the calculation within QCDSR on the light cone in Ref. 84] give 
go- Diri'mV) — 12.5 ± 1 and (7_B*S7r(0) = 29± 3 ^. The latter are small compared to lattice data or the experimental 
determination of go* Dni'mV) by the CLEO Collaboration. 
From our results we obtain the ratio 

R = 9B^BA0) fB^VM^ ^ ^ Q QQj. (54^ 

go'DTviO) Jd' \/Mb 

in good agreement with HQS that predicts a value of 1 with I/toq corrections appearing in next to leading order j6(l | 
^. Our result in Ea. (|54|l is also in agreement with the one obtained combining lattice data for Jb' and Jd- from 



Corresponding to the emission of a tt"*" 

Values for both coupUng constants obtained pri or t o 1995 within different approaches can be found in Ref. 1841 and references therein 
Note the strong couphng constant used in Ref. l60l is given in terms of ours as {qh'Ht: /vr)/(2Mjf» ) with H = B,D 
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Ref. |5J|, for gs'S-wi'^) from Ref. jsj, and the experimental CLEO Collaboration data for go* D-ni^"^) from Ref. |80l |. 
In this case one gets 

^ls.p.-La«. = 1-26 ±0.36 (55) 
where we have added errors in quadratures. A calculation using light cone QCDSR gives 

R\qcdsr = 0.92 (56) 



VI. CONCLUDING REMARKS 



Our analysis of leptonic decay constants of mesons with a heavy c or 6 quark shows that in a nonrelativistic 
calculation the equality /yMy = fpMp is satisfied within 2%. This equality is expected in HQS in the limit where 
the heavy quark masses go to infinity. The nonrelativistic result suggests that the HQS infinite mass limit sets in 
already at the rric scale, something that is not supported by lattice data for D mesons where one finds deviations as 
large as 20% from the above equality. 

One also finds problems in the semileptonic B ^ D and B — > D* decays. We have seen how the h^{w) form 
factor of the B ^ D decay and the hA2{w) form factor of the B D* decay are not reliably calculated in the 
nonrelativistic quark model where one finds large deviations from the HQET relations of Ea. l|22|l . We have tried 
to remedy this failure by evaluating the Isgur-Wise function from a form factor whose calculation we trusted in the 
quark model, (w) for the B ^ D decay, and /i^^ (w) for the B D* decay. Those Isgur-Wise function were later 
used together with HQET constraints in Ea. H22|l to recalculate all other form factors. The two Isgur-Wise functions 
thus determined show an overall reasonable agreement with lattice data but in both cases the slope seems to be 
too small. This deficiency multiply its effects when one goes to larger w values and as a consequence the quantities 
F]j{w) \Vcb\ and F^t^w) \Vcb\ go above experimental data for w values larger than 1.2, and for any reasonable value 
of \Vcb\- A failure of some kind is expected in a nonrelativistic calculation as w increases. For w = 1.2 the three 
momenta of the final meson amounts to 66% ot its mass and relativistic corrections in the wave function could start 
to be important. On the other hand we believe our results are sound at zero recoil [w — 1). That enables us to 
obtain \Vcb\ = 0.040 ± 0.006 from the B ^ D decay, in perfect agreement with the value \Vcb\ = 0.040 ± 0.005 
previously obtained by us from the analysis of the Ab — s- Ac semileptonic decay, and in good agreement with a recent 
determination \Vcb\ = 0.0414 ± 0.0012 ± 0.0021 ± 0.0018 by the DELPHI Collaboration. The experimental situation 
concerning the B — *■ D* reaction is not so clear as different experiments give values for Fd- (w) \ Vcb\ which are hardly 
compatible. From DELPHI Collaboration data we would get \Vcb\ = 0.040± 0.003 in agreement with the above result. 

Finally we have made use of PCAC to evaluate the strong coupling constants gB'Sir and go'DTr- Our results are 
larger than experimental data or the results provided by lattice and QCDSR calculations. In this case the final meson 
is nearly at rest and one would expect a nonrelativistic calculation to perform better. The main difference with the 
other observables analyzed is that here the two active quarks are the light ones. The discrepancies might hint at a 
possible sizeable renormalization of the axial coupling for light constituents quarks. On the other hand the value for 
the ratio R in Eq. (|54|l agrees with the HQS prediction and with the one evaluated using a combination of lattice and 
experimental data, being also close to a QCDSR determination. 
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APPENDIX A: MATRIX ELEMENTS FOR THE LEPTONIC DECAY OF PSEUDO SCALARS AND 

VECTOR MESONS 



The matrix element needed for the evaluation of the pseudoscalar decay constant is given by 



4^0^^ (0) 



p,o) 

/ NR 



Sl ,S2 



\st,h),{s2,f2) 



(P) 



(-1)^ 



{27r)-2^2EfAp)2EfJp) 
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1 - 



(p ) + m/ J (%,(p)+m/J 
AEfAp)Eh{p) 
\P? \ 



where we have used the fact that Vsj{p) = "f^U-sjip)- 
Similarly for the vector meson case 



4V^(o) 



K 00 



NR 



{2w)^^2Ef,{p)2Ef,ip) 

Vs2j2{p)l3Us^,fA-p) 

(-1)^-^^ 

f'^''''i2n)i^2EfAp)2EfM 

U-.2j2ip)l3l5UsiJii-p) 



= v-3fd^Y.€'kis2j2^(p)- 

Sl,S2 



4Ef,ip)Ef,ip) 

1^12 



1 + 



HEf, (p) + ruf, ) {Ef, (p) + m/J ) 



APPENDIX B: EXPRESSION FOR THE V^^Qql) MATRIX ELEMENT 

The expression for y^(|g |) is given by 



= ^2MB2En{\q\) ^jD,-\q\k 



B, 



NR 



\l2E,{\q^ + p)2Ei,{p) 



nic + mf^ 



where /2 represents a light u or d quark. 
Going a little further we have 



V%\q\) = ^2MB2ED{\q\) J d'p ^ (<^ifj^ 



(Al) 



(A2) 



(Bl) 



Ec{\q\k + p) + mA [Eb{p)+mb 



I 



^E^{\q\k+p) Eb{p) 



1 + 



\p?+pAq\ 



y \Ec{\q\k+p) + mc\ [Eb{p)+mb\ j 



(B2) 



In the case of equal masses rrib = rric and for |^| = {iv = 1) we will obtain 



(B3) 
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Similarly 



V^{\q\) = -V2MB2ED{\q\) / d'p — 



liD) 



\ 



47r {J^'f^^^ mc + mf, 
Ec{\q\k + p) + mA (Eb{p)+mb 



Pz + \q\ 



Pz 



4Ee{\q\k+p) Ebip) 



Ebip) + mb 



Ec{\q\k+p) +mc 



(B4) 



APPENDIX C: EXPRESSIONS FOR THE k'*' AND A[*^ MATRIX ELEMENTS 



The expressions for the V^*^ and matrix elements are given by 



(*) 



(|9 I) = ^2MB2En^{\q\) {D*,X- \q\ k 



B, 



NR 



2^e(|9|fc+p) 2^6(p) 



"1 

Msi,c(-|5| fc-p)7/iW4,6(-p) 



= V2MB2ED>m ^Jd*,X -\q\k 



NR 



B, 



NR 



2£e(|«|^ + p) 2i;6(p) 



^ Usi,c{-\q\ k-p) 7,^75 Us'^,b{-P) 



So that 



\ 



+ 19*1 



AE,{\q\^ + p) Ebip) \Ebip)+mb E,{\q\i + p)+m, 



(CI) 



(C2) 



■^^2MB2En^{\q\) d'p^i C/.'(l ^ 



\ 



Ee{\q\k + p) + rric^ (^Eb{p) + rub^ 
AE^{\q\k+p) Eb{p) 



I 



1 + 



2pl-\p\^-Pz\q\ 



(C3) 



= -i^2MB2Ej,,m J d'p i- ^^7m^^ |g|^ + p| )) ^ C.dPD 
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Ec{,\q\k + p) + rricj yEb{p) + mt 
4E^{\q\k+p) Eb{p) 



Pz + \q\ 



Eb{p) + mb E^^^^\% + p) + 



rUc 



(C4) 



■i^2MB2En*{\q\) j d^p i- (^c^^^^^ 



\ 



Ec{\q\k + p) + mcj yEb{p)+mb 
4Ec{\q\k+p) Eb{p) 



rric + mf. 
I 



\q\k^p\)] <^if)/l-- 



1 + 



1pl-\p\^ ^p,\q\ 



Ec(\q\k-\-p)-\-mA \ Ebi$)^mb 



(C5) 
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